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ABSTRACT 

The probability that an exoplanet transits its host star is high for planets in close orbits, but drops off rapidly 
for increasing semimajor axes. This makes transit surveys for planets with large semimajor axes orbiting bright 
stars impractical, since one would need to continuously observe hundreds of stars that are spread out over the 
entire sky. One way to make such a survey tractable is to constrain the inclination of the stellar rotation axes 
in advance, and thereby enhance the transit probabilities. We derive transit probabilities for stars with stellar 
inclination constraints, considering a reasonable range of planetary system inclinations. We find that stellar 
inclination constraints can improve the transit probability by almost an order of magnitude for habitable-zone 
planets. When applied to an ensemble of stars, such constraints dramatically lower the number of stars that 
need to be observed in a targeted transit survey. We also consider multiplanet systems where only one planet 
has an identified transit, and derive the transit probabilities for the second planet assuming a range of mutual 
planetary inclinations. 

Subject headings: methods: analytical — planetary systems — stars: oscillations (including pulsations) — 
stars: rotation — surveys 



1. INTRODUCTION 

Transiting exoplanets around the brightest stars in the sky 
will be a major provider of science in the coming years. As the 
focus of the field moves toward the characterization of tran- 
siting planets, exoplanets orbiting bright stars will provide the 
best opportunities for precise follow-up observations. This is 
particularly true for any Earth-analogs discovered around the 
brightest Sun-like stars. The Kepler mission is surveying rel- 
atively dim stars (V = 9 to V = 16) for Earth-analogs, but 
these will be difficult to observe from the ground. Detect- 
ing a transiting Earth-analog around a bright (V < 7) Sun- 
like star would enable a multitude of observations, from the 
ground and space, that would otherwise not be possible. Tran- 
siting planets around the brightest Sun-like stars are therefore 
of prime scientific interest. 

The brightest Sun-like stars have, nevertheless, yet to 
be comprehensively surveyed for transiting planets, Earth- 
analogs or otherwise. The reason is that they are spread out 
over the entire sky, which makes a practical survey of these 
stars very difficult. Since they are so widely distributed, any 
survey of the brightest Sun-like stars would have to be a tar- 
geted survey that examined one star at a time. Such a sur- 
vey wo uld be similar in concept to the M Earth survey of M- 
dwarfs ( Nutzm an & Charbonne au 2008), the N2K radial ve- 
locity survey (Fi scher et alj 120 05), or the TERMS targeted 
transit survey dKane et al.l 120091) . Unfortunately, given the 
hundreds of stars that would need to be surveyed, the aver- 
age transits probabilities, and the typical window functions 
for long period Earth-analogs, a targeted survey searching for 
long period transiting exoplanets would require either hun- 
dreds of telescopes or hundreds of years of observing time. 

In the course of considering a new space-based transit 
search for Earth-analogs, this problem motivated us to con- 
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sider ways in which a targeted survey of the brightest stars 
could be feasibly executed. One solution is to determine the 
inclination of the rotation axes of the target stars. From an- 
gular momentum considerations, we would expect any orbit- 
ing planets to lie close to the stellar equator. By only target- 
ing those stars with stellar inclinations near 90° to our line of 
sight, the number of stars that need to be surveyed could be 
lowered dramatically. 

We have therefore calculated the effect of stellar inclination 
constraints on transit probabilities. We account for a range 
of possible planetary system inclinations with respect to the 
stellar equator, and apply these probabilities to an ensemble 
of stars. We also examine the number of stars that need to 
be observed to reach various statistical confidence levels of 
detection. For reasonable assumptions, this allows for a 85% 
reduction in the number of stars that would need to be ob- 
served in a targeted survey. Furthermore, we consider the 
case of multiplanet systems where only one planet has a de- 
tected transit, and derive the transit probability of the second 
planet for various spreads of the mutual inclination angle. We 
conclude by discussing ways in which to measure stellar in- 
clinations, followed by the outline of a potential space-based 
survey. 

2. TRANSIT PROBABILITIES 

2.1. Background 

We begin by reviewing the transit probability for a sin- 
gle star under the assumption that the planet's orbital in- 
clination is randomly and evenly distributed over all possi- 
ble orientations. We will assume for now that the plane- 
tary orbit is coplanar with the stellar equator. In this case, 
the transit probability for p lanets in circular orbits is R*/a 
(Borucki & Summers] 1 19841) . the ratio of the stellar radius to 
the semimajor axis of the planetary orbit. The a priori transit 
probability of R*/a may be derived by considering the angu- 
lar momentum vector of a planetary orbit in particular orbital 
orientation, the complete set of which describe a sphere in 
space. The probability of seeing any particular set of orien- 
tations is the fractional area on the sphere that encompasses 
that set's angular momentum vectors. Again, if we assume for 
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FIG. 1. — Diagram of the relevant angles, ip is the angle that the stellar 
equator makes with the plane of the sky, A is the angle the planetary orbit 
makes with the stellar equator, and £ is the observed orbital inclination of the 
planetary orbit from Earth. 

the moment that the planetary orbit is coplanar with the stellar 
equator, the orbital inclination i is the same as the stellar incli- 
nation ip (Figure 1). The probability of star having a particular 
value for its stellar inclination then follows the distribution 



fq,(ip) = sin(V0. 



(1) 



Transits will occur if the planetary orbit has an inclination 
between (90° - 0) < i < (90° + ff), where the angle 6 is equal 
to 

6 = arcsin (~~^J > @) 

and is the maximum orbital inclination for which a planet will 
show a transit. Note that we have assumed the planet radius 
to be much smaller than the stellar radius. 

The transit probability for randomly distributed circular or- 
bits is the fraction of orbits that lie within this orbital inclina- 
tion range, and is given by 



J 90 o_ e sin(V0 dip 



f™°° sin(V>) dip 



= cos(90°- 



R* 



(3) 



which is the result quoted earlier. Note that in these and the 
following calculations we will continue to assume that the ra- 
dius of the planet is much smaller than the radius of the star. 

2.2. Transit Probabilities with Stellar Inclination Constraints 

When there is information about the orientation of the plan- 
etary system through measurement or indirect assumptions, 
the transit probability can be calculated b y treating the a bove 
distrib ution for ip as a Bayesian prior dBayesI 117631: ICoxl 
1946). This allows us efficiently to slot in stellar inclination 
measurements to our knowledge of a system's orientation. In- 
deed, the calculation of transit probabilities is almost a perfect 
subject for Bayesian techniques, since we have a rigorously 
defined prior distribution — a luxury that many problems do 
not have. 



Assuming that the planetary orbit is coplanar to the stellar 
equator, the probability distribution arising from a measured 
stellar inclination angle, ip m , can be treated as a conditional 
distribution. The resulting posterior distribution will be 



/*(V#m) = 



sin(ip)f(ip m \ip) 
foo°° sin(ip)f(ip m \ip) dip 



(4) 



The transit probability will be the integral of the posterior dis- 
tribution over the range of transiting stellar inclinations, 



P V = h(iP\iP m )dip. 
J90°-e 



(5) 



As an example, take a Gaussian measurement of the form 
ip m ± a. This implies a conditional distribution of 



fc-iusMm\lp) [ 

and a posterior of 



/Gauss(V#m) = ■ 
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(7) 



If we have made an inclination measurement of ip m = 90° ± 5° 
with Gaussian uncertainties, then for an Earth-like planet with 
R*/a = 1 /215 the transit probability will be 4.25%. This is 9 
times greater than the transit probability before the inclination 
measurement was made (0.47%). For hot Jupiters with R*/a = 
1/10, the transit probability is enhanced to 74.8% from 10%. 

Alternatively, consider a box-like, uniform, distribution of 
the measured angle. This could occur if we had a "binary" 
measurement that only revealed if the stellar inclination was 
above or below a certain angle. In this case, the conditional 
distribution from the measurement of ip m is the combination 
of two Heaviside step functions 

fboMmH)=^-(H[ip + (-ipm + <j)] H[-ip + (ip m +a)]), (8) 
La 

The posterior will be 

/box(V#m) = 

^(H[iP + (-ip m + a)] H[-iP + (iP m +a)]) 



2a 



2a 



KH[ip + (-ip m +a)]H[-iP + (.iP m + a)])diP 



(9) 



For a similar measurement of tp m = 90° ±5°, the transit prob- 
ability for an Earth-like planet at R*/a = 1/215 is 5.3%, an 
increase of a factor of 1 1 .4 over the transit probability before 
making the measurement, hot Jupiters, on the other hand, are 
close enough to the star that they will always transit across a 
star with a measured stellar inclination of ip m = 90° ±5°. 

2.3. Transit Probabilities for a Range of Planetary System 
Inclinations 

We now consider the added complexity of inclination in the 
planetary system itself. In the preceding discussion we as- 
sumed that any orbiting planets are coplanar with the stellar 
equator. From the solar system we know that this is not neces- 
sarily the case. The ecliptic is inclined by 7.155° to the Sun's 
equator, and the rest of the solar system's planets are scattered 
within several degrees of the Ecliptic. We will define the in- 
clination of the planetary system to the stellar equator — or 
any other reference plane — as the variable A. The angle the 
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stellar equator makes with the plane of the sky we denote as 
the stellar inclination ip, and we will to refer to the observed 
angle the planetary orbit makes with the plane of the sky as 
the orbital inclination i (Figure 1). 

Factoring in the planetary inclination of the exoplanetary 
system, the orbital inclination (i.e., as seen from Earth) will 
bei = ip — \, and its distribution will be governed by the joint 
distribution of i/j and A. Since ip and A are independent vari- 
ables, their joint distribution is simply the product of their 
individual distributions: fj = fy ■ /a, where fj is the joint dis- 
tribution, is the stellar inclination distribution, and f\ is 
the planetary inclination distribution. We will substitute ip - i 
for A in the distribution f\ to help the calculations. Making 
the substitution, the distribution of the angle i is then 



/i(0= / h(ip\ipm)fA(ip-i)dip. 



(10) 



The function //(/) is called the marginal probability distri- 
bution for the orbital inclination. To calculate the transit prob- 
ability of a system, //(/) becomes the integrand in Equation 



r 90° +8 , ..... 

V_e //(') di 
J o fiW di 



(11) 



Here P tr is the probability that a single star will show a tran- 
sit, assuming that there is a planet in orbit. Equations ( TTOb 
and ( fTTT i may also be used on multiple planet systems: if one 
planet is observed to transit the parent star, its orbital inclina- 
tion may be used as a reference plane in place of the stellar in- 
clination ip m , and an assumed mutual inclination distribution 
may be substituted for the planetary system inclination A. The 
transit probability of a second planet in the system may then 
be calculated. The single-star transit probability may also be 
applied to an ensemble of stars with random stellar inclina- 
tions. 

2.4. Eccentric Planetary Orbits 

We have assumed up to this point that any potential planets 
are on circular orbits. To calculate the transit probability for 
eccentric planetary orbits, consider the opening angle 8 which 
describes the maximum angular difference between our line 
of sight and the observed orbital inclination, i, that allows a 
planetary transit. In the case of a circular orbit 9 is defined as 
in Equation (f2]). In the more general case, 8 is 



' = arcsin 



A' 



(12) 



where r is the planet-star separation at the time of transit. For 
an eccentric Keplerian orbit, this implies that 



: arcsin 



R* 1 + esin(w) 
a 1 — e 2 



(13) 



Here e is the orbital eccentricity and lo is t he argument of 
periastron. This h as been noted previously bv lBarnei d2007l) 
and lBurkd d2008l) . 

In the context of transit probabilities with stellar inclination 
constraints, the effect of an eccentric planetary orbit will be 
to change the integration limits in Equation dlT~b . The exact 
effect will depend upon the values of e and lo, or more gener- 
ally upon the distribution of e and lo in exoplanetary systems. 
In the situations we are interested in, when there is no radial 
velocity evidence for a planet, we expect the distribution for 




FIG. 2. — Eccentricity effects on the transit probability for an Earth-analog, 
as a function of the measured stellar inclination angle. Note that these are 
values averaged over an assumed uniform distribution for the argument of 
periastron. We have used a measurement uncertainty of 5° and planetary 
inclinations spread uniformly within 7.5° of the stellar equator. 

the argument of periastron, to, to be evenly distributed. This 
implies that on average any orbital eccentricity will act as a 
"boost factor" to the angle 8 of 1/(1 —e 1 ). 

Figure 2 demonstrates how this boost factor affects the 
transit probability of an Earth-analog, as a function of the 
measured stellar inclination angle. We have averaged over 
an assumed uniform distribution for the argument of perias- 
tron, and used a measurement uncertainty of 5°. The planets 
are also assumed to have planetary inclinations evenly spread 
within 7.5° of the stellar equator. At e = 0.3, the boost factor 
to 9 (and, to a very near approximation, the transit probability) 
is 1.10. An eccentricity of 0.5 gives a boost of 1.33. 

Note that for hot Jupiters and Earth-analogs — the two 
cases we treat subsequently — we expect the planets to be in 
nearly circular orbits, with e < 0.1. An eccentricity of e = 0.1 
provides a boost factor of 1.01 to the angle 9. Fifty out of the 
57 transiting hot Jupiters with measured eccentricities have 
e < 0.1, and although the eccentricity distribution of Earth- 
analogs is poorly constrained and outside the scope of this 
paper, we note that Earth has an eccentricity of e = 0.017. We 
therefore only consider circular orbits in Section 3.2, when we 
calculate the benefits stellar inclination constraints provide in 
transit searches for hot Jupiters and Earth-analogs. 

2.5. Transit Probabilities for an Ensemble of Stars 

One of our goals is to determine the number of stars one 
must survey to find a transiting planet, given a fixed uncer- 
tainty in the measured stellar inclinations and an adopted 
range of planetary inclinations. To do this, we first need to 
determine the size of our initial target list, «i, for stellar in- 
clination measurements. This will set the probability P ntr that 
there is the desired number of transiting planets in the stars 
that we are considering. For instance, if we measured the stel- 
lar inclination of just 10 stars there would be a low probability 
that a transiting Earth-analog is in our sample, but measuring 
the stellar inclination of several thousand stars would make it 
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near certain. 

Second, we also want to calculate the fraction of the initial 
target stars that will actually need to be observed photometri- 
cally for transits. This fraction will set the odds that we will 
observe one of the transiting planets hidden within our origi- 
nal set of stars, which we denote as P b s . Since we are truly 
interested in only those stars with stellar inclinations near 90°, 
the size of this fraction will depend directly on the precision of 
the stellar inclination measurements. It will also be sensitive 
to assumptions made about the distribution of A, the plane- 
tary system inclination. Later, it will be useful to delineate 
the fraction of the initial stars that we will observe photomet- 
rically by the angle <j>. We then will observe all of the stars 
with measured stellar inclinations within cf> degrees of 90°. 

We first turn to describing P m , the probability that there is 
the desired number of transiting planets within our initial n 
target stars. We will calculate f ntl - using the binomial distribu- 
tion. Previous studies have instead used the expectation value 
to determine the number of stars that one needs for a transit 
detection: if hot Jupiters have a 10% transit probability, then 
a survey will need to look at ten stars to statistically expect a 
single detection. We use the binomial distribution to perform 
a similar calculation. Unlike using the expectation value, the 
binomial distribution allows us to explicitly place a probabil- 
ity that the survey will detect the desired number of transiting 
planets. This calculation will be based upon the a priori tran- 
sit probability denoted as P lr fi = R*/a, since we are trying to 
determine how many stars we need to start with in our initial 
sample to yield enough transiting planets. We will denote the 
binomial distribution for s successes in n independent trails, 
each with success probability p, as the function 



Bi(n,s,p)- 



(pYH-pT 



(14) 



If we start with n\ initial target stars, consider the probabil- 
ity ^det.i that there is exactly one transiting planet in the initial 
stars, and that we detect it. This probability will be the prod- 
uct of the odds P ntr i that there is exactly one transiting planet 
in the n\ initial stars, and the odds that we will observe it pho- 
tometrically, P b s . Using the binomial distribution to calculate 
the chance that there is exactly one transiting planet in our 
initial stars, we can write 

— * ntr,l ^obs = Bi(ni,l,P tr .o)P b s . (15) 

Now let us take the case of exactly three transiting planets 
among our initial stars. What is the probability that we will 
detect at least one of them? Again we may use the binomial 
distribution to calculate the probability P ntr 3 that there are ex- 
actly three transiting planets, but now we will also use a sum 
of the binomial distribution to determine the probability that 
we observe at least one of these planets. We take the sum be- 
cause we will regard observing one, two, or all three of the 
planets as successfully fulfilling our criteria that we observe 
at least one planet. We therefore add the probabilities for all 
three of these scenarios, and 

3 

^det,3=Bi(« 1 ,3,n-.o) ^Bi(3,£,P obs ). (16) 

k=\ 

The general case, and the one that we are ultimately inter- 
ested in, is the one wherein we wish to know the probability 
that we will detect at least n tl transiting planets, but do not 
specify the exact number of transiting planets in our initial n; 
stars. To calculate this overall probability we must sum the 



individual probabilities that we will detect at least n a plan- 
ets given an exact number of transiting planets in our ti\ stars. 
This similar to how we summed the probability of observation 
in Equation ( TToT l. In terms of our initial target stars, the proba- 
bility of photometrically observing one transiting planet, and 
the number of transiting planets that we wish to detect, the 
probability that we will detect at least that many planets is 
therefore 

■Pdet(«i 1 ^Obs 1 W tl ) = ]T Bi(«„;,P tr , ) ^Bi(j,fe,P obs ) . 

j=n u \ &="lr / 

(17) 

We are still left with the calculation of P b s , the indepen- 
dent probability that we will observe each of the transiting 
planets within our initial «; stars. We will only photometri- 
cally observe a fraction of the initial ti\ target stars for tran- 
sits. These will be the stars with measured inclinations within 
<f> degrees of 90°, and the size of the fraction will be deter- 
mined by the desired level of P b s . A survey will only pho- 
tometrically observe stars with measured stellar inclinations 
between 90° - 4> < ip m < 90° + 4>. Now consider that within 
the total number of stars, that will have their stellar in- 
clinations measured we expect that the measured inclinations 
will be randomly distributed around the sky, and so the distri- 
bution of tp m will go as sm(ip m ). The expected number of stars 
that will lie within <fi degrees of 90° will therefore be 



J 9 oo_^, sin(V> m ) d% 
/o^sin^m) dip m 



«obs = «i 



= msin(0). (18) 



To find the number of stars, n„f, s that will need to be observed 
for transits, we will need to calculate n\ and <f>. Note that 
we are considering the general case and therefore leave out 
survey-specific considerations that should be included for a 
more accurate yield estimate. For example, we will assume 
that every star has a planet of the type we are looking for, 
whereas the true frequency will depend upon the period or 
mass range that the survey is most sensitive to. We also do 
not consider the effect of non-central transits or variations 
in stellar-type, since these will both be directly tied into the 
signal-to-noise and magnitude limits of a given survey. For a 
more detailed discussion on how these factors can affect the 
yields of transit surveys, see|Beattv & Gaudij d2008l) . 

To calculate cj), the range of measured stellar inclinations 
that will be considered for transit observations, first consider 
the transit probability as a function of the measured stellar in- 
clination of the target star. This will allow us to calculate the 
probable distribution of the transiting planets as a function of 
the measured stellar inclination, and thereby figure what stel- 
lar inclination range will have to observed to recover a given 
fraction of the transiting planets within our sample of «; stars. 

Including the inclination of the planetary systems as per 
Equation ( TTOb . the transit probability will be 



Ar(V'm) = 



J90°-e Jo° 



fcaus S (i>\i>m)fA(^-i) dip di 



Jo° 8 ° Jo° 8 ° fGa U s S {^Hm)fh{lp-i) dtp di 



(19) 



where fGauss(ip\ipm) is the posterior distribution of the stellar 
inclination. We now multiply P tr by the probability that mea- 
suring a star yields a given stellar inclination angle, sin(V> m ). 
Doing so will give us the probability that measuring the stellar 
inclination of a star yields a certain value, and that with a mea- 
sured stellar inclination ip m the star will then show a transiting 
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FIG. 3. — Diagram of a notional multiplanet system displaying the angles 
referenced in the text. The angles % and i c are the orbital inclinations of the 
planets as observed from Earth. A(, c is the mutual inclination angle between 
the two orbits. 

planet, 

P* m (1>m) = sin(V> m )f'tr(V'm)- (20) 

The distribution P^ m is, physically, the expected number of 
planets — normalized to unity — that we will expect to see 
transiting stars with a given measured stellar inclination ip m . 

The fraction of transits that we will recover is dependent 
upon how much of this distribution is observed for transits. If 
we observe all of the stars with measured stellar inclinations 
within an angle <f> of 90°, then 



^obs 



Jg^ sin(V>m) Ptr(tpm) dlp m 
J I80 sin(V' m ) Pn(lpm) d^ m 



(21) 



For a given value of P \, s , we may solve for the angle cf> that 
defines the outer limit of the subsample observed for transits. 

We can therefore calculate n\ and <f> for a given survey by 
specifying the desired number of transiting planet detections, 
and the desired values of P ^. This allows us to determine 
«obs, the number of planets that will have to be observed for 
transits, as well as the probability that we will detect the de- 
sired number of transiting planets, 



fdet(«i,-Pobs,"tr) = ^ I Bi ("i,;\-Ptr,()) ^ Bj(j, fc, P obs ) 
j=n a \ k=nu 



(22) 



3. RESULTS 



3.1. Probability for a Second Planet to Transit 

One direct application of our probability calculations is in 
determining the transit probability of a second planet in a sys- 
tem in which one planet is already observed to transit the par- 
ent star. From angular momentum considerations, we have 
been using the stellar equatorial plane as the reference plane 
for our assumed distributions for the inclination of the plan- 
etary systems, /a(A). In the case of a multiple planet sys- 
tem we may instead use the plane of the transiting planet's 




a priori probability - 



A 



be, max 



FIG. 4. — Transit probability of HAT-P-13c, as a function of the maximum 
mutual inclination of the two planets (A(, c max ). This assumes that the orbit of 
planet 'c' is evenly distributed within A(, c max degrees of planet'b's orbit. For 
reference, the a priori transit probability of R t /r c is also plotted. 

orbit as the reference. Consider the planetary system shown 
in Figure 3. We may measure the orbital inclination 4 m of 
the transiting inner planet through the transit photometry. As- 
suming Gaussian uncertainties, this will give us a distribution 
// b ('b|'b,m) for the orbital inclination of the inner planet. The 
outer planet will be inclined at some angle Ab c to the inner 
planet, and will have an orbital inclination with respect to the 
sky of z' c = /b + Abc ■ The mutual inclination of the two planets 
Abc will in turn be drawn from the distribution ,/A hl .(Ab c ). 

Inserting these distributions into Equation ( [Tol l gives the 
marginal distribution of the orbital inclination for the outer 
planet (here we replace ip m with the angle i'b as the relevant 
reference plane for the system) 



flc('c)= / // b 0'b|i6,m)/A bc 0'c-l'b) di b . 



The transit probability will be 



P* = 



J9O°-0 C JlA l c) ai c 
f 180= 

J o fl c ( l c)dl c 



(23) 



(24) 



as per Equation (fTTT i. The angle 8 C is determined by the orbital 
separation of planet 'c' and the host star at the time of their 
conjunction, C = arcsinCR^/rc). 

As an example, consider the HAT-P-13 system (Bakos et al. 
120091) . As of this writing 4 only the inner planet, HAT-P-13b, 
has been observed to transit. What is the probability that the 
outer planet. HAT-P-13c, also transits? iBakos etal] d2009t) 
measure the orbital inclination of the inner planet to be «b ; m = 
83.4° ± 0.6°. At the time of its conjunction with HAT-P-13, 
the outer planet is at a distance of r c /R* = 82.1 ±6.1 stellar 
radii from the star. This allows us to determine // b (i b |'b.m) 
and 9 C for the HAT-P-13 system. For the mutual inclination 
of the two planets, we assumed that the orbital inclination of 

4 2009 September 
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planet 'c' was evenly distributed within ±Ab c .max degrees of 
the orbital inclination iy, of planet 'b'. Figure 4 shows the 
transit probability of HAT-P-13c as a function of Abc.max- 

The probability that planet 'c' transits is very dependent 
upon what we assume is a reasonable range of mutual incli- 
nation in the HAT-P-13 system. For reference, all of the Solar 
System planets are within 3.4° of the Earth's orbit — except 
for Mercury (at 7°). If the mutual inclination of the two plan- 
ets orbiting HAT-P-13 is within 3.4°, then the outer planet 
will not transit. If the two planets are misaligned by up to 
7°, then the transit probability for planet 'c' is 7%. The max- 
imum transit probability of 8.5% occurs if we assume that 
the two planets may be inclined within 8° of each other. As 
the assumed spread in mutual inclination increases, the transit 
probability will fall back to the a priori value of 1.2%. 

The transit probability for HAT-P-13c is therefore at most 
8.5%. 

3.2. Benefits for Target Selection 

We now demonstrate how using stellar inclination measure- 
ments and the enhanced transit probabilities (Section 2) can 
aid in the target selection of transit surveys. As illustrative 
cases, we will calculate how many stars need to be observed 
in a survey looking for hot Jupiters, and for a separate survey 
searching for planets within the habitable-zone. In these ex- 
amples, we will make the simplifying assumption that every 
star has either a hot Jupiter or habitable-zone planet in orbit, 
at distances of R*/a = 1/10 or R*/a = 1/215, respectively. 

We assume that the inclinations of the planetary systems are 
distributed in two ways. For the hot Jupiter s, we use the plane- 
tary in clination distribution determined by Fabrycky & Wind 
(2009) from an ensemble of 11 Rossiter — McLaughlin mea- 
surements of spin — orbit alignment. The authors found that 
aside from the XO-3 system 5 the hot Jupiters they consid- 
ered had planetary inclinations distributed according to a 
Rayleigh distribution with a width parameter of 6.6°. Exo- 
planets within the habitable-zone may not follow this same 
planetary inclination distribution. We use a uniform distribu- 
tion of planetary inclination within 7.5° of the stellar equator; 
Earth has an planetary inclination of 7.155° to the Sun's equa- 
tor. Including either spread of planetary inclinations into the 
calculations of transit probabilities acts to spread out the prob- 
ability of transit, and make stars with stellar inclinations far 
from 90° more likely to show transits. At the same time, the 
spread of planetary inclinations makes stars with measured 
stellar inclinations near 90° less likely to show transits. As- 
suming a measurement of tp m = 90° ± 5°, and that planetary 
inclinations are uniformly distributed within 7.5° of the stel- 
lar equator, then the transit probability for a habitable-zone 
planet at a distance of R*/a = 1/215 star drops from 4.25% to 
3.08% as compared to assuming the orbit is coplanar with the 
stellar equator. Conversely, a star with measured stellar incli- 
nation of ip m = 80° ± 5° has its transit probability increased 
from 0.59% to 1.11% 

As our first example, take a survey for hot Jupiters around 
solar-type stars. We will require a 95% probability that the 
survey detects at least a single transiting planet. From Equa- 
tion ([T7| > this means that Pdet = 0.95. Although we may ar- 
bitrarily set P b s and «i, as shown in the appendix the time 
required to complete a hot Jupiter survey is minimized if we 
set P b s = 0.9814. To find the number of stars needed in the 

5 Which appears to be highly misaligned, possibly because of XO-3b's 
migration history. 



initial sample we must then solve 

P det = 0.95 = jr UiimJ,^-) ]TBi(M,0.9814)j . (25) 
j=\ \ k=i ) 

We must therefore have w; = 29 stars in our initial sample. 

We next want to know how many stars out of these 29 will 
actually have to be observed photometrically. That is, how 
many of the initial targets with measured inclinations near 90° 
will we need to look at for transits? We will assume that all 
of the stellar inclination measurements have Gaussian uncer- 
tainties of 5°. The transit probability for a hot Jupiter can be 
calculated for various orientation measurements of the form 
ip m ±5° by using Equation ( TTOb to account for the inclination 
of the planetary system: 



Pti(4>m) = 



fgo°-5j3° fo°° fc ai i SS (^\^m)fA(tp-i) dip di 



fo° Jo° 8 ° fGauss(^\lpm)fK(i>-i) dip di 



(26) 



hot Jupiters at a distance of R*/a = 1/10 will show transits up 
to a maximum angle of 9 = 5.73°. The angle <p that defines 
our observed subsample solves 



JofW, sin(V>m) Ptrilpm) dlp m 

P obs = 0.9814= ^ 

Jo sin(V> m ) P a (ipm) dtp m 



(27) 



and is = 24.01° 

To detect at least one hot Jupiter, we must therefore photo- 
metrically observe 



n obs = 29 sin(24.01°)= 11.80 



(28) 



stars that will have measured stellar inclinations within 24.01 
of 90°. This will give us a probability of Pa et = 95% of detect- 
ing at least one hot Jupiter. The top panel of Figure 5 shows 
how the number of stars that need to be observed varies as a 
function of the stellar inclination measurement precision for 
various confidence levels. The top panel also shows how the 
fraction of the initial target list that will need to be observed 
varies with measurement precision. The lower limits in both 
cases set b y the spread in the distr ibution of /a (A) as deter- 
mined by EEyc^&Wni dHH). 

In a notional survey for habitable-zone planets, we will also 
require a 95% chance of detecting at least one transit. The cal- 
culations are similar to those for the survey for hot Jupiters, 
except that we will assume that f\ (A), the distribution of plan- 
etary system inclinations, is evenly distributed within 7.5° 
of the stellar equator. Otherwise, we keep the requirement 
that fdet = 0.95, and optimize P Q b s in the same way. The time 
needed to complete a survey for Earth analogs is minimized 
if we set P b s = 0.7746, which means the initial sample size 
will be «; = 830 stars. Proceeding in the same way as the 
hot Jupiter survey, we will need to observe stars with mea- 
sured stellar inclinations out to <p> = 8.12° — corresponding 
to n bs = 1 17.2 stars that will need to be observed photomet- 
rically. The lower two panels of Figure 5 show how n bs 
and how the fraction of initial stars that need to be observed 
change as a function of measurement precision, assuming the 
planetary inclinations are within 7.5° of the stellar equator. 

Unlike the case of the hot Jupiters, there is no statistically 
compelling information regarding the inclination distribution, 
/a, of Earth-analogs. Since the fraction of the stars that would 
need to be observed photometrically is tied to this distribu- 
tion, Figure 6 illustrates the effect of changing the maximum 
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FIG. 5. — Number of stars that would need to be photometrically observed 
in a survey for different detection confidence levels. Note that this assumes 
every star has a corresponding planet. 

spread of the planetary inclination on the curves from the bot- 
tom panel of Figure 5. In the above calculations we used a 
uniform distribution for f\ out to a maximum of 7.5°, but 
here we also plot the fraction of stars that will need to ob- 
served against stellar inclination measurement precision for 
a maximum of 15° and 0° (i.e. perfectly aligned). Note 
the difference between the three curves becomes most pro- 
nounced with higher precision (lower uncertainty) inclination 
measurements as the underlying distribution for f\ becomes 
more dominant. Under our assumption that the uncertainty in 
ip m is ±5°, the fraction of stars that need to be observed is 
0.1 1, 0.14 and 0.21 for maximum inclinations of 0°, 7.5° and 
15°, respectively. 

In both of our example surveys, the number of stars that 
need to be observed photometrically is dramatically lower 
than the initial sample size. The exact difference will depend 
upon the precision of the stellar inclination measurements; we 
have assumed Gaussian uncertainties of 5°. The uncertainties 
themselves depend upon the method used to measure stellar 
inclinations. The two most prominent ways to measure stellar 
inclinations are either through spectroscopic measurements of 
stellar vsin(?/>) or through asteroseismic observations. 

4. DISCUSSION 

4. 1 . Spectroscopic v sin(^) Inclination Measurements 

Spectroscopic measurements of vsin(?/>) allow us to mea- 
sure stellar inclinations indirectly. If the true rotational period 
of the star can be identified through photometric variation — 
or other means — and we have an estimate for the radius 
of the st ar, we will be able to ca lculate the val u e of s injijj). 
Recently IWinn et alJ (I2007I) and lArentoft et all (I2008I) have 
measured sin(^) for two different stars. The Arentoft group's 
work is also particularly illustrative of the potential pitfalls 
in any attempt to measure sm(ip): they observed photomet- 
ric variation on Procyon with a period of 10.3 ± 0.5 days, but 
found it much more likely that the true rotational period of 
Procyon is twice that value. 
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FIG. 6. — Fraction of stars from the initial target list that need to observed 
photometrically as a function of the uncertainty in measuring the stellar in- 
clination. The three curves plotted correspond to differing maximums for the 
spread of the planetary inclination distribution. 

The use of vsinQip) measurements to constrain stellar in- 
clinations is limited by the uncertainties in the measure- 
ment of vsin(^), the stellar radius, and the stellar rotational 
period. Spectroscopic measurements of vsin(^) for solar- 
like stars typically h a ve fractional uncerta inties of 15-25% 
( Kepp ens et ail [19951 iTerndrup et all 120021) . The precision 
may be increased by using extremely high resolution spec- 
tra (R sal00,000; Carney et al. 2008) to disentangle the line 
variations caused by rotation and turbulence. To measure 
the rotational period of a star, one would either need clearly 
identifiable photometric variation, or would have to observe 
a star for sev eral months to identi fy the rotation period in a 
periodogram (Meibom et al. 2009). Together with fractional 
uncertainties in the stellar radius and rotational period, it is 
therefore time intensive to measure sin(?/») to better than 15- 
30%. 

Using measurements of sin(^) to constrain stellar inclina- 
tions is complicated by the flatness of the sine function near 
90°. Consider that arcsin(0.9) = 64° and that arcsin(0.995) = 
84°. Any determination of sin(^) will therefore need to be 
extremely precise to yield usable constraints on angles near 
ip = 90°. As noted previously, this level of precision would 
require specialized, time-intensive observations. However, 
measurements of vsin^) can also be used to identify and 
eliminate from consideration stars with stellar inclinations far 
from 90°. Assuming a 15% fractional error on sin(^), one 
would be able to identify 28% of stars as unsuitable. 

4.2. Asteroseismic Inclination Measurements 

Stellar inclinations may also be measured directly through 
precise asteroseismological measurements of solar-like oscil- 
lations. These 5 minute acoustic oscillations in the stellar 
photosphere can be described as spherical harmonics with 
the harmonic numbers n, /, and m. On a non-rotating star 
the 21 + 1 ;7?-modes are degenerate and lay on top of one 
another in frequency-space. As the angular velocity of the 
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star increases, the m-modes un dergo rotationa l — splitting, and 
pull apart from one an other (lLedouxlll951l) . Additionally, 
iGizon & Solankil d2003l) show that the relative power in each 
of the m-modes will depend upon the stellar inclination. By 
measuring the magnitude of the splitting and the relative 
power of the split modes, it is possible to determine the an- 
gular velocity and stellar inclination of a star that undergoes 
solar-like oscillations. 

Asteroseismic observations c an be conducted using e ither 
photometry or spectroscopy. IGizon & Solankil d2003l) and 
iBallot et"al] d2008l) describe the theoretical basis and the ex- 
pected uncertainties in measuring stellar inclinations using 
photometric asteroseismological observations. Problemati- 
cally for surveys of Sun-like stars, both papers calculate that 
stars with angular velocities near solar are extremely chal- 
lenging targets. In Gizon & Solankj] d2003l) the formal error 
on measurements of the stellar inclination can be very large as 
a star's angular velocity approaches that of the Sun. Neverthe- 
less, photometric asteroseismology is being conducted from 
space by the CoRoT, Kepler, and MOST missions, and on the 
ground by numerous observers. To date, the only photometric 
detection of rotational splitting in a star other than the Sun has 
been accomplished photometri cally using the CoRoT space- 
craft {Appourchaux et al. 2008); though the authors note in- 
consistencies in their data as compar ed to earlier observations 
of the same star (M osser et all 2005). 

Spectroscopic asteroseismological observing programs 
have also multiplied. While none have identi fied rotational 
splittin g in an unevolved main sequence star, iBouchy et alJ 
(2005) observed rotational splitting in the acoustic spectrum 
of the G3IV-V star p Arae, and sever al other groups hav e 
come close to an identification (see, e.g jBazot et al. (2007)). 
One of the most ambitious sp ectroscopic observing collab- 
orations is the SONG project (Grunda hl et al.1 120 08). which 
aims to build several dedicated telescopes spaced in longitude 
around the world to allow for continuous observing. Spec- 
troscopic observations should theoretically provide more pre- 
cise stellar inclinations measurements than photometry, since 
the widths and shapes of the absorption lines provide addi- 
tional information about the photosphere that is not present 
in luminosity variations. Work is currently underway to char- 
acterize the exact stellar inclination precision that can be ex- 
pected from spectroscopic asteroseismology (T. Campante & 
H. Kjeldsen, 2009, private communication). 

4.3. Practical Applications 

Precise measurements of stellar inclinations are a step to- 
ward a practical transit survey of the brightest Sun-like stars 
in the sky. One can envision a space-based survey for transit- 
ing exoplanets targeted at these stars. The idea of manufactur- 
ing and launching a suite of nanosatellites, each with a single 
telesc ope and targe ted at an individual star is currently under 
study (Seager 2008). The goal of the study would be to design 
and build a space telescope that fits within a lOx 10x30 cm 3 
triple CubeSat spacecraft, and to take advantage of the grow- 
ing number of piggyback launch opportunities to place these 
telescopes into low-Earth-orbit. The present major challenge 
is the required pointing stability to achieve a high enough pho- 
tometric precision to detect the transit of an Earth-analog on 
a low-mass satellite (< 5kg). 

The duration of the survey would be largely set by the num- 
ber of stars that would need to be observed photometrically 
for transits. Each star would need to be observed for at least a 
year to cover the full orbital period of an Earth-analog; though 



it may be possible to assign individual spacecraft multiple 
targets. For the initial stellar inclination measurements, us- 
ing spectroscopic vsin(V>) inclination determinations would 
require photometric observations spread over at least a stel- 
lar rotation period — on the order of a month for solar-type 
stars. Asteroseismic inclinations would be, on average, faster 
to come by: SONG estimates that their network, with six ded- 
icated observing stations across the globe, could measure the 
inclination of up to a few dozen stars over one year (H. Kjeld- 
sen, 2009, private communication). 

In terms of actual telescope time, the photometry for 
vsin(-0) inclination determinations would require on average 
10 minutes a night over the course of a month, for a total 
of 5 telescope-hours per star over that month. The astero- 
seismic inclination measurements, using SONG's estimates, 
would require about 10 telescope-days per star. Our fiducial 
Earth analog search would require inclination measurements 
on 830 stars, which corresponds to about 170 telescope-days 
for vsin(^) measurements, and about 8400 days for a 6 tele- 
scope SONG network. The photometric transit survey por- 
tion would take approximately 120 telescope-years. Note that 
even with the time required for the inclination measurements, 
this is one seventh the time required for a comparable blind 
survey. 

The idea of a targeted space-based transit search is made 
feasible by the enhanced transit probabilities that are achiev- 
able using stellar inclination constraints. One possible mis- 
sion concept would be to use ground-based vsin^) measure- 
ments to eliminate a third of the possible targets from con- 
sideration. Asteroseismology conducted from the ground and 
from orbiting triple CubeSat spacecraft could then be used to 
assemble a final target list for photometric observations. De- 
pending upon the precision of the asteroseismology, the de- 
sired number of detections and the desired confidence level of 
achieving this many detections, this would reduce the number 
of stars that would need to be observed from a few thousand to 
a few hundred. A related proposal, usi ng only ground -based 
asteroseismology has been described by Beatty (2009). 

5. SUMMARY 

Transit surveys of the brightest solar-like stars are made dif- 
ficult by the distribution of these targets across the entire sky. 
This makes a traditional point-and-stare photometric survey 
unworkable, since each target star would need to be observed 
individually. Such a survey of the brightest (V < 7) stars 
would require either a prohibitive number of telescopes or a 
prohibitive amount of observing time, especially if it is tar- 
geting Earth-analogs in stellar habitable zones. This problem 
motivated us to consider the effect of stellar inclination mea- 
surements on transit probabilities and the number of stars that 
need to be observed in a photometric survey to statistically 
expect a detection. 

We derived the transit probability for stars, individually and 
in an ensemble, with constraints on the stellar inclinations and 
with assumptions about the range of planetary system incli- 
nations. This derivation involved several steps. First, the stel- 
lar inclination constraint was treated as the conditional dis- 
tribution in Bayes theorem and combined with our prior as- 
sumption of randomly oriented orbits. Second, the planetary 
inclination was included via a joint probability distribution 
together with the stellar inclination constraints. These steps 
completed the probability distribution for a single star. These 
calculations also allowed us to compute the transit probability 
for the second planet in a multiplanet system when the other 
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planet has already been observed to transit. In this case, we 
used the orbital plane of the transiting planet in place of the 
stellar inclination as the reference plane for the system, and 
the mutual inclination of the two planets in place of the plane- 
tary system inclination. An ensemble of stars was treated us- 
ing the binomial distribution and the distribution of expected 
transit detections to determine the number of stars required to 
be observed, and the probability that doing so would yield a 
detection. 

In summary, our first result is the transit probability of a 
single planetary system that has measured stellar inclination 
constraints. Assuming inclination measurements accurate to 
5°, we find that the transit probabilities for typical hot Jupiters 
may be increased to 74.8%, and the transit probabilities for 
Earth-analogs may be increased to 4.25%. These calculations 
may also be applied directly to finding the transit probability 
of the second planet in a multiplanet system where one planet 
has already been observed to transit. For the specific case of 
HAT-P-13, we find that the transit probability for the outer 



planet is between 0-8.5%, dependent on what is the assumed 
spread of the mutual inclination. Our second result is the esti- 
mated number of stars needed to be observed for transits given 
a specified required number of planet detections and desired 
probability of achieving those detections. Assuming 5° un- 
certainty on the stellar inclination measurements, we would 
need to look at 120 stars to have a 95% chance of detecting 
more than one Earth-analog. We would have a 50% chance 
of detecting more than 3.8. This is one seventh the number of 
stars that a blind transit survey would need to look at for the 
same yield. 
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and discussion. This work was supported in part by the NASA 
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APPENDIX 
OPTIMIZING Pobs AND M 

The overarching goal in optimizing the P b s and n\ is to minimize the length of time needed to complete the survey. Recall that 
we start with the initial set of n\ target stars whose stellar inclinations we measure. The exact value of n\ will determine, through 
the statistics of the binomial distribution, the probability P m that a given number of transiting planets are within our initial target 
stars. We then observe a fraction of these initial stars that have stellar inclinations close to 90° photometrically for transits. The 
size of this fraction sets the probability P Q b s that we observe photometrically the transiting planets hidden within our initial stars. 
As described in Equation d2lT >, the size of this fraction is dependent upon the precision of the stellar inclination measurements; 
the more precise the measurements, the smaller the fraction may be for a fixed P b s . 

Consider our notional survey to discover transiting Earth-analogs. As described in Section 4.3, it would take on average a month 
per star to measure stellar inclinations — roughly enough time to see one complete rotation of a Sun-like star. The photometric 
transit observations would take at least one year per star. Assuming that every star has an Earth-analog, that stellar inclinations 
may be measured to ±5° and that planetary inclinations are evenly distributed within 7.5° of the stellar equator, Figure 7 shows 
the effect of varying P„b s and n\ on the time needed to complete the survey — while keeping the overall detection probability 
constant at 0.95. As the number of initial targets increases, the corresponding fraction (and hence P b s ) of those stars that need to 
be photometrically observed decreases. The minimum amount of time needed to complete the survey occurs when n; = 830 and 
P obs = 0.7746. 

While we do not show the optimization of a hot Jupiter survey in Figure 7, the method is exactly that same as for an Earth 
analog survey. The only major difference is that the photometric observations take about 10 days to complete. Since the time 
needed to measure the stellar inclination remains one month, this pushes a hot Jupiter survey to minimize the number of initial 
targets that must be observed for inclination measurements. If we wish a 95% chance of detecting at least one hot Jupiter, then 
this means we will have to set n; = 29 and P b s = 0.9814 — assuming that every star is orbited by a hot Jupiter. 
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FIG. 7. — Optimization of the P a j, s and m for an Earth-analog search that will have a 95% confidence of at least one detection. The minimum amount of time 
required to complete the survey (normalized to unity in the figure) occurs when P a \, s = 0.775 and n\ = 830. 



